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Rotational transitions in molecular hydrogen collisions are computed. The two most recently 
developed potential energy surfaces for the H2— H2 system are used from the following works: 
1) A.I. Boothroyd, P.G. Martin, W.J. Keogh, M.J. Peterson, J. Chem. Phys., 116 (2002) 666, 
and 2) P. Diep, J.K. Johnson, J. Chem. Phys., 113 (2000) 3480; ibid. 112, 4465. Cross 
sections for rotational transitions 00^20, 22, 40, 42, 44 and corresponding rate coefficients are 
calculated using a quantum-mechanical approach. Results are compared for a wide range of 
kinetic temperatures 300 K < T < 3000 K. 

1. INTRODUCTION 

The interaction and collision properties of hydrogen molecules, and hydrogen molecular iso- 
topes has been of great theoretical and experimental interest for many years [EEHH1I11E1IE1I3IE1 
II3[ini[IIl[Il[I3[IS[Iia[ia[I3 Because of the low number of electrons 

in the H2— H2 system this is one of the few four-center systems for which the potential energy 
surface (PES) can be developed with very high precision. Therefore H2+H2 is a benchmark col- 
lision used for testing several dynamic methods such as: semiclassical [El; quantum-mechanical 
[ I19| or wave packet ['20 studies. The system may also be useful in improving our understanding 
of fundamental processes in few-body molecular dynamics. Additionally, the H2+H2 elastic and 
inelastic collisions are of interest in combustion, spacecraft modeling and in clean renewable 
energy. Hydrogen gas, in particular, has generated much interest as a energy supplier, see for 
example [EE]- 

The hydrogen molecule plays an important role in many areas of astrophysics [ 261 1271 12*51 129j . 
It is the simplest and most abundant molecule in the universe especially in giant molecular 
clouds. Energy transfer involving H2 molecules governs the evolution of shock fronts [ 1261 127] and 
photodissociation regions (PDRs) in the interstellar medium. Collision-induced energy transfer 
between H2 molecules and between H2 and other atoms/molecules is related to an important 
astrophysical processes, which is the cooling of primordial gas and shock wave-induced heating in 
the interstellar media. To accurately model the thermal balance and kinetics of such important 
systems one needs accurate state-to-state rate constants. 

It is well known, that experimental measurements of quantum state resolved cross sections 
and rates is a very difficult technical problem. On the other hand accurate theoretical data 
requires precise PESs and reliable dynamical treatment of the collision processes. The first 
attempt to construct a realistic full-dimensional ab initio PES for the H2— H2 system was done 
in works [ 1101 lllj . and the potential was widely used in the framework of variaty of methods 
and computation techniques. 
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Because of the immense theoretical and practical benefits associated with the recent hydrogen 
fuel issues, the H2— H2 system has been reinvestigated and an accurate interaction potential from 
the first principles has been developed in work [ 30 . However, in this work the Diep and Johnson 
potential energy surface (DJ PES) was extrapolated only for the rigid rotor monomer model 
of H2— H2. On the other hand two extensive studies of the H2— H2 PES have been reported 
by Boothroyd et al., In these studies the potential energies have been represented at 

6101 and 48180 geometries respectively with a large basis set at the multireference configuration 
interaction level. The earlier 6101 points were fitted to a six dimensional many-body expansion 
form in work [I14j. 

In this work we present a comparative study of the global BMKP and DJ PESs for collisions 
of rotationally excited H2 molecules. The scattering cross sections and their corresponding rate 
coefficients are calculated using a non reactive quantum-mechanical close-coupling approach. In 
the next section we will shortly outline the method. Our results and discussion are presented 
in Section 3. Conclusions are provided in Section 4. Atomic units (e=m e =7i=l) are used 
throughout this work. 

2. METHOD 

In this section we will briefly present the close-coupling quantum-mechanical approach we 
used to calculate our results, specifically the cross sections and rates in collision of a hydrogen 
molecule with another hydrogen molecule. The Schrodinger equation for a ab + cd collision in 
the center of a mass frame, where ab and cd are linear rigid rotors is 

where is the momentum operator of the kinetic energy of the collision, R is the collision 
coordinate, Myi is a reduced mass of the pair of two-atomic molecules (rigid rotors in this 
model) ab and cd: M\2 = {ma + m^^c + m^/irna + m^ + mcA-md), ^1(2) are reduced masses 
of the targets: 1^1(2) = f^a{c) m b{d) I '( m a(c) + m b(d)), ^1(2) are the angles of orientation of rotors 
ab and cd, respectively, J is total angular momentum quantum number of abed system and M 
is its projection onto the space fixed z axis, V(fi,r2, R) is the potential energy surface for the 
four atomic system abed. 

The eigenfunctions of the operators in (^Q) are simple spherical harmonics Yj iUli (f). To 

solve the equation Q the following expansion is used [E] 



TjJM ( r>\ 

Wur 2 ,R)= E j R & 2 L^,r-.R). (2) 

JMjxj 2 ji2L 

where channel expansion functions are 

^^12^1^2, R) = E ^jimii2m2^/i2mi2«m^imi(^l)^2m 2 (^2)^Lm(-R), (3) 
ra\m2m\2Vti 

here j\ + ji = ii2> 3X2 + L = J, mi, m,2, m\2 and m are projections of ji, j\2 and L 
respectively. 

Substitution of (J2J) into Q provides a set of coupled second order differential equations for 
the unknown radial functions U^ M (R) 

( d 2 L(L + 1) 



dR 2 R 2 



+ k 2 a U J a M {R) = 2M 12 E / < <f>i M (h,r 2 , R) 



a' 

|V(fi,r 2 ,i?)|^/ M (fi,f 2! i?) > UiF (R)dhdr 2 dR, (4) 
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where a = (jij2ji2^)- We apply the hybrid modified log-derivative- Airy propagator in the 
general purpose scattering program MOLSCAT [ 32 to solve the coupled radial equations (J3J). 
Additionally, we have tested other propagator schemes included in MOLSCAT. Our calculations 
revealed that other propagators can also produce quite stable results. 

The log-derivative matrix is propagated to large i?-intermolecular distances, since all exper- 
imentally observable quantum information about the collision is contained in the asymptotic 
behaviour of functions U£ M (R — > 00). The numerical results are matched to the known asymp- 
totic solution to derive the physical scattering S"-matrix 

Ui ~ 5 aa ,e-^^-^/ 2 ))-f^V /2 ^ a ,e-^^-aV2)) j (5) 

R-^+co \k aa i ) 

where k aa > = 2M\2(E + E a — E a i) 1 / 2 is the channel wavenumber, E a ^ a i^ are rotational channel 
energies and E is the total energy in the abed system. The method was used for each partial 
wave until a converged cross section was obtained. It was verified that the results are converged 
with respect to the number of partial waves as well as the matching radius, R ma x, for all channels 
included in our calculations. 

Cross sections for rotational excitation and relaxation phenomena can be obtained directly 
from the 5-matrix. In particular the cross sections for excitation from jxj 2 — » jij 2 summed over 
the final m^m^ and averaged over the initial m\m2 are given by 

{2 3l + l){2 3 2 + l)k aa , Jmj , 2LL , 

S J {j' 1 ,j2,j'x2L';h,32,h2,L-E)\ 2 . (6) 

The kinetic energy is e = E — £?ui(ji + 1) — £? 2 j 2 (i2 + 1)> where -61(2) are the rotation constants 
of rigid rotors ab and cd respectively. 

The relationship between the rate coefficient fe^-^' j' (T) and the corresponding cross section 
a j 1 j 2 ^j' 1 j' 2 (Ekin) can be obtained through the following weighted average 

knJ^(T) = ^(^f jT (7) 

where e = Et t a i—Ej 1 —Ej 2 is precollisional translational energy at the translational temperature 
T and e s is the minimum kinetic energy for the levels j\ and j 2 to become accessible. 

3. RESULTS 

In this section we present our results for rotational transitions in collisions between par a /para- 
hydrogen molecules: 

H 2 (ii = 0) + H 2 (j 2 = 0) - H 2 (j0 + H 2 (j 2 ). (8) 

We apply the newest PESs from the works [I3U] and [ I31j . The first one, DJ PES, is constructed 
for the vibrationally averaged rigid monomer model of the H 2 — H2 system to the complete basis 
set limit using coupled-cluster theory with single, double and triple excitations. A four term 
spherical harmonics expansion model was chosen to fit the surface. It was demonstrated, that the 
calculated PES can reproduce the quadrupole moment to within 0.58 % and the experimental 
well depth to within 1 %. 

The bond length was fixed at 1.449 a.u. or 0.7668 A. DJ PES is defined by the center-of-mass 
intermolecular distance, R, and three angles: 6\ and 62 are the plane angles and (pi2 is the 
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relative torsional angle. The angular increment for each of the three angles defining the relative 
orientation of the dimers was chosen to be 30°. There are 37 unique configurations for each 
radial separation when the symmetry of the H2— H2 system is considered. In previous works 
calculating the potential a much smaller set of angular configurations designed to represent the 



represents an expansion on Legendre polynomials [[3j. 

The second potential, BMKP PES, is a global six-dimensional potential energy surface for two 
hydrogen molecules. It was especially constructed to represent the whole interaction region of 
the chemical reaction dynamics of the four-atomic system and to provide an accurate as possible 
the van der Waals well. The ground state and a few excited-state energies were calculated. The 
new potential fits the van der Waals well to an accuracy within about 5% and has an rms error 
of 1.43 millihartree relative to the 48180 ab initio energies. For the 39064 ab initio energies that 
lie below twice the H2 dissociation energy BMKP PES has an rms error of 0.95 millihartree. 
These rms errors are comparable to the estimated error in the ab initio energies themselves. In 
the six-dimensional conformation space of the four atomic system the conical intersection forms 
a complicated three-dimensional hypersurface. The authors of the work [|^ mapped out a large 
portion of the locus of this conical intersection. 

The BMKP PES uses cartesian coordinates to compute distances between four atoms. We 
have devised some fortran code, which converts spherical coordinates used in Sec. 2 to the 
corresponding cartesian coordinates and computes the distances between the four atoms. In all 
our calculations with this potential the bond length was fixed at 1.449 a.u. or 0.7668 A as m 



Now we will present our results for the elastic and inelastic integral cross sections and rate 
coefficients for the collision (JHJ). As far as astrophysical applications are concerned, we are 
particularly interested in the pure rotational transitions of the H2 molecules. 

A large number of test calculations have been done to secure the convergence of the results 
with respect to all parameters that enter into the propagation of the Schrodinger equation 
This includes the intermolecular distance R, the total angular momentum J of the four atomic 
system, Ni v i the number of rotational levels to be included in the close coupling expansion and 
others (see the MOLSCAT manual [I32|1. 

We reached convergence for the integral cross sections, a(Eki n ), in all considered collisions. 
In the case of DJ PES the propagation has been done from 2 A to 10 A, since this potential is 
defined only for the specific distances. For the BMKP PES we used r m j n = 1 A to r max = 30 
A. We also applied a few different propagators included in the MOLSCAT program. 

Table 1 represents the convergence test results with respect to J max-, the maximum value of 
the total angular momentum, for both the BMKP and DJ PESs. The calculations are limited 
to just three values of energy, for the simpler basis set: Jij2=00, 20 and 22, and from lowest 
to highest within the considered range of energies. As can be seen the results are stable for 
the range of kinetics energies, when J max is increased from 80 to 90. In all our subsequent 
production calculations we use J max = 80. 

It is important to point out here, that for comparison purposes we don't include the compen- 
sating factor of 2 mentioned in [ 9 j. However, in Fig. 2 and in our subsequent calculations of 
the rate coefficients, kjj'(T), the factor is included. 

In Table 2 we include the results of our test calculations for the various rotational levels j'u'2 
included in the close coupling expansion. In these test calculations we used two basis sets: 
j lj2 =00, 20, 22, 40, 42 with total basis set size N M = 13 and jij 2 =00, 20, 22, 40, 42, 44, 60, 




DJ PES. 
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Figure 1. Rotational state resolved integral cross sections for j\ = 32 = — > j[ = 2,j' 2 = 0, 
ji =32 = 0^ j{ = 2,f 2 = 2, ji = j 2 = -» = 4,j 2 = 0, ji = j 2 = -» = 4, j 2 = 2 and 
Jl = J2 = — ¥ Ji = 4, j 2 = 4 calculated with the BMKP and DJ PESs (the compensating factor 
of 2 is not included). 



62 with Ni v i = 28. One can see that the results are quite stable for the 00— >20 and 00^22 
transitions and somewhat stable for the highly excited 00^40 transition. Nontheless, for our 
production calculations we used the first basis set. 

The objective of this work is to make reliable quantum-mechanical calculations for different 
transitions in P-H2+P-H2 collisions and provide a comparative study of the two PESs. The 
energy dependence of the state-resolved integral cross sections o 'j l j 2 -»j' 1 j' 2 {Ekin) for the j\ = 32 = 
— > j[ = 2,j' 2 = and j\ = J2 = — > j[ = 2,j' 2 = 2 rotational transitions are represented 
in Fig. 1 (upper plots) for both the BMKP and DJ PESs respectively. These channels for the 
most part influence the total cross-section of the H2+H2 collision. As can be seen both PESs 
provide the same type of the behaviour in the cross section. These results are in basic agreement 
with the recent calculations of work [EOj, where the BMKP PES was also applied, but using a 
time-dependent quantum-mechanical approach. Our calculation show, that DJ PES generates 
higher values for the cross sections, by up to 50%. 

The integral cross sections for inelastic collisions: j% = 32 = — * j[ = 4,j 2 = 0, j\ = ]2 = 
-> j[ = 4, j 2 = 2, and jt = j 2 = j{ = 4, j' 2 = 4 for BMKP and DJ PESs are presented in 
Fig. 1 (bottom plots). The cross sections are very small at energies less than 0.25 eV, since at 
lower kinetic energies these transitions are closed by the corresponding energy barriers. 

We would like to point out here, that the j\ = 32 = — > j[ = 4, j' 2 = 2 cross section becomes 
larger than j% = 32 = — > j[ = 4:,j' 2 = if the collision energy is greater than ~0.6 eV. In this 
energy range it is hence more likely that the second diatom is also excited when the first diatom 
makes the — > 4 transition. However the j\ = 0, 32 = — > j[ = 4, j 2 = 4 cross section is very 
small over the entire energy range considered. The DJ potential energy surface again provides 
higher results in the corresponding cross sections. One can note, that for both potentials the 



Table 1: Convergence of the total cross sections (10 -16 cm 2 ) for transitions 00— >20 and 00— >22 with respect to the maximum value of the 
total angular momentum J max in the H2— H2 system. Here a B and a D are the cross sections calculated with BMKP and Diep and 
Johnson [|3U] PESs respectively. 
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Table 2: Convergence of the total cross sections (10 _16 cm 2 ) for transitions from 00— >20, 22, 40 with respect to the number Ni v \ of the 
levels to be included in the basis set of the H2— H2 system. Here a B and a D are the cross sections calculated with BMKP [|^ and Diep 
and Johnson [ 30 PESs respectively. 

20, 22, 40, 42 (N M = 13) Basis set: jij 2 =00, 20, 22, 40, 42, 44, 60, 62 (N M = 28) 
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Relative kinetic velocity [m/s] 



Figure 2. Integral elastic cross sections calculated with the BMKP and DJ potentials. The 
experimental measurements are those of Bauer and co-workers [E] (the compensating factor of 
2 is included). 




Temperature [KJ 



Figure 3. Temperature dependence of the state-resolved thermal rate constants for the j\ = 

h = -> fi = 2,32 = M = 2 >i2 = 2,fi = ±,32 = and 3i = 4,3*2 = 2- The results for the 
DJ and BMKP PESs are given in solid and dashed lines, respectively. The diamonds are the 
theoretical data of work [llfij. 
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Table 3 

Rate coefficients fcoo^jj'(^) (cm 3 s~ 1 ) calculated with the DJ PES for rotational transitions 
00 — > 02,22 and 40 in comparison with other theoretical [E] an d experimental [|23] data. 
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rotational inelasticity is dominated by the 00 — > 20 transition. 

The integral cross sections for an elastic P-H2+P-H2 collision computed with both PESs are 
depicted in Fig. 2. The two cross sections are in reasonable agreement over a wide range of 
energies. Again the BMKP PES generates higher values. The difference becomes even larger 
at lower kinetic energies. In the figure we also provide experimental data from work The 
theoretical and the experimental results are in reasonable agreement with each other. This fact 
indicates, that the spherical parts of BMKP and DJ PESs are close in shape, which is a very 
important attribute for PESs. 

The differences in the cross sections between the two potentials are also reflected in the state- 
resolved thermal rate constants, as shown in Fig. 3. Again, the BMKP PES underestimates the 
rate constant for the 00— * 20 transition, and overestimates those transitions of higher rotational 
levels. The near perfect agreement for the 00— ► 40 transition is likely accidental. Because of the 
highly averaged nature of the rate constant, the difference is not as conspicuous as in the cross 
sections. We have also plotted in the same figure the rate constants reported previously in work 
[ITS'], where Schwenke's PES [ HO] was used. 

Finally, our rate constants for 00 — > 20, 22, 40 rotational transitions calculated with the DJ 
PES at lower temperatures are listed in Table 3 together with the theoretical calculations of 
Flower et al. [E] an d recent experimental results from work [US]. As can be seen our rate 
constants /coo^2o(^) an d £;oo-+22(r) are close to those of work [Ej. The experimental results 
are higher by about 60% for 50 K and 30 % for 110 K. 

4. SUMMARY and CONCLUSIONS 

A systematical study of the state-resolved rotational excitation cross sections and rates in 
molecular par a- /par a- hydrogen collisions is completed. A test of convergence and the results 
for cross sections and rate coefficients using two different potential energy surfaces for the H2— H2 
system have been obtained for a wide range of kinetic energies. 

Although our calculations revealed, that both PESs can provide the same type of behaviour 
in regard to cross sections and rates, there are still significant differences. The DJ potential 
overestimates by about 20-40 % the results at even relatively larger kinetic energies. This is 
especially true in regard to 00— >20 rotational transition, where significant differences at around 
300 K are seen in Fig. 3. 
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Considering the results of these calculations one can conclude that subsequent work is needed 
to further improve the H2— H2 PES. Detailed calculations including rotational- vibrational basis 
set and comparative analyses using both potentials at low and very low kinetic energies for 
0-H2/0-H2 and P-H2/0-H2 excitation-deexcitation collision processes currently are in progress in 
our group. 
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